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The current status of the experimental measurements and theoretical predic- 
tions of the anomalous magnetic moment of the muon is briefly reviewed. 
The emphasis is put on the evaluation of the hadronic contribution to as 
it has the largest uncertainty among all Standard Model contributions. The 
precision of the hadronic contribution is driven by the input e~^e~ data pre- 
dominantly from the vr+vr" channel. Including the latest experimental data on 
e^e~ annihilation into hadrons from CMD2 and SND for the tt~^tt~ channel and 
BaBar for multihadron final states, the updated Standard Model prediction 
disagrees with the measurement dominated by BNL by 3.3 standard deviations, 
with the theoretical precision exceeding the experimental one. 



1 Introduction 

For a charged elementary particle with 1/2 intrinsic spin such as muon, its magnetic dipole 
moment /2 is aligned with its spin s as: 



where q = ±e is the charge of the particle in unit of the electron charge and g is the gyromagnetic 
ratio. In the classic Dirac theory, g = 2. In the Standard Model (SM), quantum loop effects 
induce a small correction, which is quantified by = {g^ — 2)/2, the so-called anomalous 
magnetic moment or the magnetic anomaly. 

There has been a long history in measuring and calculating a^. In particular the steadily 
improving precision of both the measurements and the predictions of and the disagreement 
observed between the two have led the study of one of the most active research fields in 
particle physics in recent years. 

The paper is organized as follows. In section [2l the measurement history and the current 
world average value of are presented. In section [3l different components of the SM contri- 
butions to are reviewed. Section [H is reserved for discussions followed by conclusion and 
prospects in section [5l 
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2 The Measurement of 

A compilation of the major experimental efforts in measuring over the last five decades 
is given in Table [1] (a modified version of Table 1 from a recent review article Starting 
from the experiment at the Columbia-Nevis cyclotron, where the spin rotation of a muon in a 
magnetic field was observed for the first time, the experimental precision of has seen constant 
improvement first through three experiments at CERN in the sixties and seventies and more 
recently with E821 at the Brookhaven National Laboratory (BNL). The current world average 
value reaches a relative precision of 0.54 ppm (parts per million). 



Table 1: Measurements of the muon magnetic anomaly a^, where the value in parentheses stands for 
either the total experimental error or the statistical and systematic ones. 



Experiment 


Beam 


Measurement 




Columbia-Nevis(1957)EI 
Columbia-Nevis(1959)'^ 


^+ 


g = 2.00 ±0.10 

0.001 i3;^jj2) 


12.4% 


CERN 1(1961)^ 
CERN 1(1962)13 
CERN 2(1968)E1 
CERN 3(1975)E1 
CERN 3(1979)121 


^+ 
/^^ 


0.001 145(22) 
0.001 162(5) 
0.00116616(31) 
0.001 165 895(27) 
0.001 165 911(11) 


1.9% 
0.43% 
265 ppm 
23 ppm 
7.3 ppm 


BNL E821(2000)E1 
BNL E821 (2001)1121 
BNL E821(2002)1II] 
BNL E821 (2004)1121 


^+ 


0.001 165 9191(59) 
0.001 165 920 2(16) 
0.001 165 920 3(8) 
0.001 165 9214(8)(3) 


5 ppm 
1.3 ppm 
0.7 ppm 
0.7 ppm 


World Average(2004)li^LL3j 




0.001 165 920 80(63) 


0.54 ppm 



The muon magnetic anomaly in all modern experiments is determined by the following 
method. For an ensemble of polarized muons which are moving in a storage ring in a highly 
uniform magnetic field B (perpendicular to muon spin and orbit plane) and a vertically focusing 
quadrupole field E, the frequency difference LOa between the spin procession lOs and the cyclotron 
motion uJc is described by 

^ ^ 0xE) 



UJa = UJ. 
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when B-P = E-fi = (2) 



where f3 represents the muon direction. The second term in parentheses vanishes at 7 = 29.3 
(the magic momentum) and the electrostatic focusing does not affect the spin. The key to the 
experiment is to determine frequency uJa to high precision and to measure the average magnetic 
field to equal or better precision. 

In comparison with the electron magnetic anomaly, Oe is more precisely measuredl^^ (0.7 ppb), 
but is more sensitive to new physics effects by about m'j^/m'^ ~ 40 000 because of its large 
mass value. 



3 Prediction of the Standard Model Contributions 

In the SM, the muon magnetic anomaly receives contributions from all electromagnetic 
(QED), weak and strong (hadronic) sectors and can be conveniently written as: 

^SM ^ ^QED ^ ^weak ^ ^had ^ (3) 
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Their representative diagrams are shown in Fig.[Tl which also includes two example contributions 
from new particles in supersymmetry models. Thus comparison of the precision measurement 
and theory tests the validity of the SM at its quantum loop level and probes effects of new 
physics. 




/i /i 




jj, iJ, 




Figure 1: Representative diagrams contributing to a,j. First column: lowest-order diagram (upper) 
and first order QED correction (lower); second column: lowest-order hadronic contribution (upper) 
and hadronic light-by- light scattering (lower); third column: weak interaction diagrams; last column: 
possible contributions from lowest-order Supersymmetry. 



3.1 QED and Weak Contributions 

The QED correction, which includes all photonic and leptonic (e, // and r) loops, is by far the 
dominant contribution in the SM: 



QED ^ ^0.765 857410(27) (^^^ + 24.050 509 64(43) (^^^ + 130.9916(80) (^^^ + 

+663(20) + (4) 



where the lowest-order Schwinger term (a/27r) was known since 1948 the coefficients are 
analytically known for terms up to (a/vr)^, numerically calculated for the fourth term and 
recently estimated for the fifth term^^. Using a extracted from the latest measurement^^, 
one has 

^QED ^ 718.09(0.14)5th order(0.08)5„ X 10-" . (5) 

The week contributions, involving heavy Z, or Higgs particles, are suppressed by at 
least a factor -jjr- — 4 x 10^^. At one-loop order, 



5 1/ o \2 f ml \ ( 

+ -(l-4sin^.H.) -f 



ar^'^il-loop] = — ^ 
^ ^ 8\/27r2 Lo o ^ ' yivi^ 

M2 



3 ' 3 I • - \ M^I • ' \Mjj 



(6) 



194.8 X 10-" , for sin^^ = 1 - ~ 0.223 . (7) 

z 



Two-loop corrections are relatively large and negative 

^weakp_jQQp] ^ _4o.7(i.o)(1.8) X 10-" , (8) 



3 



where the errors stem from quark triangle loops and the assumed Higgs mass range Mh = 
15014q'^ GeV. The three-loop leading logarithms are negligible, 0(10"^^), implying in total 

al^""^ = lh^{l){2) X . (9) 



3.2 Hadronic Contributions 

The hadronic contributions are associated with quark and gluon loops. They cannot be cal- 
culated from first principles because of the low energy scale involved. Fortunately, owing to 
unitarity and to the analyticity of the vacuum polarization function, the lowest-order hadronic 
vacuum polarization contribution to can be computed via the dispersion integrally using 
the ratio (s) of the bare cross section § for e'^e annihilation into hadrons to the pointlike 
muon pair cross section at center-of-mass energy y/s 

1 /a\2 K{s) 



had,LO 

3 Vvr 



dsf^R(^){s), (10) 



where K{s) is the QED kerneltIS] K{s) = (l - ^) + (1 + xf (l + ^) [ln(l + x) - x + ^ 



+ 



Inxj^, with X = jrfj^ s-^id P/j. = f 1 ^ j • The kernel function K{s) ^ ^ gives weight to 

the low energy part of the integral. About 91% of the total contribution to aJ^^^'LO jg accumulated 
at ^/s below 1.8 GeV and 73% of aj^^*^'^^ is covered by the vrvr final state, which is dominated 
by the /o(770) resonance. 

A detailed compilation of all the expe riment al data used in the evaluation of the dispersion 
integral prior to 2004 is provided in Refs.^^OlSI] gjnce then, a few precise measurements have 
been published. A list of experiments for the dominant vrvr channel is shown in Table [2j 



Table 2: A list of measurements of e'^e annihilation into hadrons in the 7r"'"7r (7) channel. 



Experiment 


^data 


Energy range (GeV) 


5(stat.) 


(5(syst.) 


DM1 (1978)1^ 


16 


0.483 - 1.096 


(6.6 - 40)% 


2.2% 


TOF (1981)123] 


4 


0.400 - 0.460 


(14 - 20)% 


5% 


OLYA (1979, 1985)1^41251 


2 + 77 


0.400 - 1.397 


(2.3 - 35)% 


4% 


CMD (1985)125] 


24 


0.360 - 0.820 


(4.1 - 10.8)% 


2% 


DM2 (1989)123 


17 


1.350 - 2.215 


(17.6 - 100)% 


12% 


CMD2 (2003)123 


43 


0.611 - 0.962 


(1.8-14.1)% 


0.6% 


KLOE (2005}ES] 


60 


0.600 - 0.970 


(0.5-2.1)% 


(1.2-3.8)% 


SND (2006)129] 


45 


0.390 - 0.970 


(0.5-2.1)% 


(1.2-3.8)% 


CMD2io„ (2006)130] 


10 


0.370 - 0.520 


(4.5 - 7)% 


0.7% 


CMD2rho (2006)131] 


29 


0.600 - 0.970 


(0.5-4.1)% 


0.8% 


CMD2high (2006)132] 


36 


0.980 - 1.380 


(4.5 - 18.4)% 


(1.2 - 4.2)% 



The vrvr data are compared in Fig. [2l Closer inspections show that the most precise measure- 
ments from the annihilation experiments SND and CMD2 at Novosibirsk are in good agreement. 
They differ however in shape with those measured by KLOE using the radiative return method 
at 

da$neI2S] 

(see Sec. 14. ip . Before this is clarified, the KLOE data are not used in some of the 
recent evaluations of al^^'^'^^. 

"The bare cross section is defined as the measured cross section corrected for initial state radiation, electron 
vertex contributio ns an d vacuum polarization effects in the photon propagator but with photon radiation in the 
final state included^I^. 
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Figure 2; Comparison of tt+tt spectral functions expressed as e'^e cross sections. The band 
corresponds to combined data used in the numerical integration. 



In addition to the dominant vrvr mode, results from the BaBar experiments are being 
produced on multihadron final states using also radiative return ESi. Benefiting from its big 
initial center-of-mass energy of 10.6 GeV, hard-radiated photon detected at large angle and high 
statistics data sample, the BaBar measurements are precise over the whole mass range. One 
example is shown in Fig. [3] in comparison with earlier measurements. 
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Figure 3: The measured cross section for e'^e — > 27r"'"27r from BaBar compared to previous measurements. 

Including these new input e'^e~ data, a preliminary update § of a^^'^'^'^ is performed and 
shown in Table El There is no new tau data since the previous evaluation, therefore the r based 
calculation is taken directly from Ref. ^^B. The evaluation using r data is made^^^l \)y relating 
the vector spectral functions from r — > 1'^-+ hadrons decays to isovector e~^e~ hadron s cross 
sections by isospin rotation. All known isospin breaking effects are then taken into accountl22l2l]_ 



''It is preliminary as some of the new e^e data used were still in their preliminary form. 
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Table 3: Summary of the • contributions from e'^e annihilation and r decays. 



Modes 


Energy [GeV] 


hoH 1 ()/ 4- \ r-i n 1(11 

Qliacl,l.U(^g+g )[10 ^'^J 


had ,{ ) / \ ri rv — 1 h 

aj^'^'^'^^(r)[10 ^"] 


Low s expansion 


2m7r 


- 0.5 


55.6 ± 0.8 ± O.lrad 


56.0 ± 1.6 ± 0.3su(2) 


-L _ 


0.5 - 


- 1.8 


449.0 ± 3.0 ± 0.9rad 


464.0 ± 3.0 ± 2.3su(2) 




2m7r 


- 1.8 


16.8 ± 1.3 ± 0.2rad 


21.4 ± 1.3 ±0.6su{2) 


27r+27r~ 


2mn 


- 1.8 


13.1 ± 0.4 ± O.Orad 


12.3 ± 1.0 ±0.4su(2) 




0.3 - 


0.81 


38.0 ± 1.0 ± 0.3rad 




0(1020) 


1.0- 


1.055 


35.7 ± 0.8 ±0.2rad 




Other excl. 


2m,r 


- 1.8 


24.3 ± 1.3 ±0.2rad 




J/V^,V(25) 


3.08- 


- 3.11 


7.4 ± 0.4 ± O.Orad 




[QCD] 


1.8- 


-3.7 


33.9±0.5qcd 




R [data] 


3.7- 


- 5.0 


7.2 ± 0.3 ± O.Orad 




R [QCD] 


5.0 - 


- oo 


9.9±0.2qcd 




sum 


2m^ 


— oo 


690.8(3.9)(1.9)rad(0.7)QCD 


710.1(5.0)(0.7)rad(2.8)su(2) 



The higher order (NLO) hadronic contributions ' involve one hadronic vacuum po- 
larization insertion with an additional loop (either photonic or leptonic or another hadronic 
vacuum polarization). They can be evaluated^S^ with the same e^e~ — > hadrons data sets used 
for a\\^'^'^. The numerical value^^Sl ^gg^j^g 

had,NLO ^ _9.79(o.09)exp(0.03)rad X IQ-'^ (11) 



where the first and second errors correspond respectively to the experimental uncertainty of the 
e'^e~ data and the radiative correction uncertainty. 

Another higher order hadronic contribution to is from the hadronic light-by-light scatter- 
ing (illustrated with the lower figure in the second column in Fig. [T]). Since it invokes a four-point 
correlation function, a dispersion relation approach using data is not possible at present. In- 
stead, calculations'^ involving pole insertions (or Goldstone boson exchanges), short distance 
quark loops and charged pion (and kaon) loops have been individually performed in a large Nc 
QCD approach. A representative value used in Ref.^^was a^^'^'^^^ = 8.6(3.5) x lO^^''. 

A new analysis^^, which taks into account the proper matching of asymptotic short-distance 
behavior of pseudoscalar and axial-vector contributions with the free quark loop behavior, leads 
to ah*^<i'LBL _ ]^3_g(|2.5) X 10~^'^. However, as pointed out in Ref.1221, several small but negative 
contributions such as charged pion loops and scalar resonances were not included in the latter 
calculation, thus in a recent update evaluation!^ of a^, the following value 

^had,LBL ^ 12.0(3.5) X 10-1° (12) 

was used|. This is consistent with the value aJ^^^diLBL _ 11(4) x 10^^^, suggested in Ref.^B. The 
uncertainty a^^'^'^^^, being the second largest one next to a)^^'^'^^, clearly needs improvement 
in the near future. 



21l36l42l43l4Hi0l 



Adding all SM contributions together, the comparison from recent evaluations 
with the measurement is shown in Fig. [H While the r data-based calculation agrees with the 
measurement within the errors, the e'^e" data-based evaluations show a deviation of around 3.3 
standard deviations. 



'^A different evaluationBS used directly the value of Ref-^SSlof 13.6(2.5) x 10" 
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Figure 4: Comparison of recent theoretical evaluations of with the BNL measurement. 



4 Discussions 



4.1 Tau Data versus e^e Data 

The r data used in the evaluation is the averaged one from the LEP experiments 
and and the CLEO experiment The data are compared in Ref.^I^ and found in 

good agreement in particular for the two most precise data from ALEPH and CLEO. These 
data are complementary as the ALEPH data are more precise below the p peak while CLEO 
has the better precision above. 

A comparison between the averaged r data and the e~^e~ for the dominant vrvr mode is shown 
in Fig. [5l The difference of 5 — 10% in the energy region of 0.65 — 1.0 GeV^ is clearly visible. 
The difference with KLOE is even more pronounced. 




0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 0.5 0.6 0.7 0.8 0.9 



s (GeV^) s (GeV^) 

Figure 5: Relative comparison of the n'^Tr~ spectral functions from e^e~-annihilation data and 
isospin-breaking-corrected r data, expressed as a ratio to the r spectral function. The shaded band 
indicates the errors of the r data. The right hand plot emphasizes the region of the p peak. 
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4.2 CVC 



An alternative way of comparing r and e^e~ data is to compare measurements of branching 
fractions B in r decays with their expectations from CVC (Conserved Vector Currect) using 
e^e~ spectral functions, duly corrected for isospin breaking effects. The advantage of a such 
comparison is that the measurements of B are more robust than the spectral functions as the 
latter ones depend on the experimental resolution and require a numerically delicate unfolding. 

The comparison for vrvr mode revealing a discrepancy of 4.5 standard deviations is shown in 
Fig. [6l Similar comparisons for decay modes t~ —>■ UT-'n'~'i'K^ and r~ — > i't-2-k~ -k^ tt^ have also 
been made and the differences with the corresponding e'^e~ data are found respectively at 0.7 
and 3.6 standard deviations. 



T decays 



Belle (preliminary) 
25.15 + 0.04±0.31 

CLEO 

25.42 + 0. 12±0.42 
ALEPH 

25.47±0.10±0.09 

DELPHI 

25.29±0.20±0.14 

OPAL 

25.44 + 0. 17±0.29 

T average 

25.40 + 0.10 



e*e CVC I Average (w/o KLOE) 

I I 24.48 + 0.18 



23.5 24 24.5 25 25.5 26 26.5 27 
B(x"^V^7t"7t°) (%) 

Figure 6: The measured branching fractions for r~ — > Vt'k'^'k^^ compared to the expectation from the 
e^e~ 7r"*'7r~ spectral function applying the isospin-breaking correction factors. 



5 Conclusion and Prospects 

The muon magnetic anomaly is one of the most precisely known quantities both experimen- 
tally and theoretically in the SM. Incorporating new e^e~ data from CMD2 and SND for tttt 
mode and from Babar for multihadronic modes, new SM determinations of have been ob- 
tained with a theoretical precision exceeding for the first time in recent years the experimental 
one. The SM prediction is found to be smaller than the measurement by about 3.3 standard 
deviations. Unfortunately one can not draw a definitive conclusion for the moment as the r 
data based prediction is in agreement with the measurement. 

Therefore it is extremely important that one clarifies the discrepancy between the e^e^ and r 
data in particular on the tttt mode. There are a number of possibilities: (1) (the normalization of) 
the e'^e data is wrong, (2) the tau data are wong, (3) both are correct but there are unaccounted 
effectsHSl 

which explain the discrepancy between the two. 
Possibility (1) may be also related to the current difference (mainly on the shape of the 
spectral functions) between CMD2/SND data and KLOE data obtained respectively from the 
beam energy scan method and the radiative return method. This difference is expected to be 
resolved soon as KLOE has more and high quality data to be analyzed. In addition, reduced 
systematics uncertainties can be achieved if the measurement is made by normalizing the tttt 
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data to instead of to luminosity using the large angle Bhabha process. 

The long awaited high precision measurement in vrvr mode from Babar using also the ra- 
diative return method will certainty help in clarifying some of the issues. 

On the tau side, further improvement on the high mass part of the spectral functions is 
expected from large statistical data samples available at B factories and a r-charm factory. 

While the leading hadronic uncertainty gets improved with the forthcoming high precision 
e~^e~ (and r) data, the next item awaiting for significant improvement concerns the uncertainty 
on the light-by- light scattering contribution (jhad.LBL 

Given the fact that the theoretical error is already smaller than the experimental one, it is 
timely to improve the latter. Indeed there is a new project BNL-E969 allowing to reduce the 
current error by more than a factor two down to 0.24 ppm. We are looking forward that the 
project gets funded very soon. 
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